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Tonks-Girardeau gas, super- Tonks-Girardeau gas, and bound 
states of one-dimensional bosons in a hard-wall trap 
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Abstract. - We investigate the Bose gas with repulsive or attractive interactions between atoms 
in the scheme of Bethe Ansatz equation in a hard wall trap. Three typical quantum phases in 
the current research of ID interacting cold atoms are clarified in terms of the energy spectrum, 
single particle density distribution and two-particle correlation function. We identify two matching 
points in the phase diagram, i.e. the TG and STG gas show the same profiles at the strongly 
interacting point — 1/7 = 0, and in the weakly interacting limit 7 = the ground states TG and 
BS join to each other smoothly. 



Recently experiments and theories have made a great 
progress in one dimensional quantum gases [T] . Thanks to 
the Feshbach resonance as well as the confinement-induced 
resonance (CIR) [2HI], the interaction strength between 
atoms can be tuned exactly. The first breakthrough in ex- 
periment is the realization of Tonks-Girardeau (TG) gases 
of bosons with strongly repulsive interaction. The 
atoms in TG gas stay in ground state and the wave func- 
tion and energy was found by Girardeau in 1960 using a 
Fermi-Bose (FB) mapping [7], leading to fermionization 
of many properties of this Bose system. The many-body 
properties were determined explicitly for a bosonic Tonks- 
Girardeau gas confined in the potential of a harmonic trap 
[8,9 or a hard wall trap [10] with a tunable 6- function bar- 
rier at the trap center and related work has been done for 
Bose-Ferimi mixture [9,11 and finite interacting strength 

nun]. 

The second important experiment is the realization of 
the super Tonks-Girardeau (STG) gas [15], where the Ce- 
sium atoms are prepared in a long-lived, strongly interact- 
ing, excited state. Unlike the strong repulsively interact- 
ing TG gas, the interaction of bosons in STG state is at- 
tractive. Although the TG gas has proved very successful 
both in theories and experiments, bosons with attractive 
interaction start to attract more and more attention only 
recently. In fact, if the interaction is finite and negative, 
the ground state is McGuire's cluster state [T61IT7] . which 
is a bound state and would decay quickly via molecular 
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Fig. 1: (Color online) The energy spectrum of N = 4 bosons 
for the ground states and one of the excited states in the full 
interacting regime. Quantum phases of TG, STG and Bound 
State are clearly classified into three branches solution of Bethe 
Ansatz Equations. 



channels. In the original proposal, the STG gas-like state 
corresponds to a highly excited Bethe state in the inte- 
grable interacting Bose gas for which the bosons acquire 
hard-core behavior [18,19 . Experimentally the STG gas 
was realized through a sudden quench of the interaction 
parameter [20]. On the other hand, the ground state of 
a strongly attractive spin-1/2 Fermi gas can be effectively 
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described by the STG gas [2"Pj . 

Exact solutions are available only in very few limiting 
cases. ID Bose gas with arbitrary repulsive interaction 
strength was solved by Lieb and Liniger |22j using the 
Bethe ansatz method in which particles are constrained 
by periodic boundary conditions on a line of length L. As 
far as the trapping of the atoms is concerned, a solution 
was given for two ^-interacting particles in a harmonic 
trap in [23] and the STG wave function for N — 2 parti- 
cles is found to be identical to the hard-sphere bose gas 
|24j given exactly in terms of a parabolic cylinder func- 
tion. It is very interesting to consider what happens to 
the above mentioned three branch solutions of ID Bose 
gas when the atoms are confined in a hard wall trap, i.e., 
under the open boundary condition. Interacting bosons 
in a hard-wall trap is one of the few integrable examples 
of interacting many-body models and its exact solution 
has been obtained with the superposition of Bethe's wave 
function for the repulsive interaction in a seminal paper 
by [25] . For a finite number of bosons and finite system 
size, boundary effects are expected to be pronounced at 
low temperature. Indeed, significantly different quantum 
effects should be exhibited by a finite number of bosons 
confined in a finite hard wall box. 

Here we study the property of a few bosons trapped in 
ID hard wall trap by numerically solving the Bethe Ansatz 
equations in the full interacting regime. We compare the 
energy of the system, the one body density matrix and 
the two-particle correlation function both in the ground 
state for both repulsive and attractive interactions and in 
the excited state for attractive interaction. The transition 
probability from the TG gas to the STG gas is examined 
through a sudden switch of interaction. Typical features 
are shown to be quite different for the three branches of 
Bethe Ansatz solution, which may serve as a guide line for 
the experiments with hard wall boundary condition. 

Consider N bosonic particles of mass m in the ID 
hard wall of width L with 8 function interaction. The 
stationary behavior of the many-body wave function 
^(xi, X2, ■ ■ ■ xn)) is described by the Schrodinger equa- 
tion 
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The interaction strength g\p> = —2h 2 /mciip, and a\u = 
a± (C — a±/a,3£)) is the ID scattering length. The tightly 
confining wave-guides in experiments are related to the 
theoretical ID models through the so-called Confinement 
Induced Resonance (CIR). Across the resonance point 
a±/a 3 £> = C = 1.0326, the g\p> can be changed from 
— oo to +oo by (i) varying the transverse width of the 
waveguide a± = y/ a / 'muj j_ where is the frequency 
of the transverse confinement, and/or (ii) tuning the 3D 
scattering length cl^d by means of a magnetically in- 
duced Feshbach resonance. In the paper we use the di- 
mensionless Tonks-Girardeau parameter 7 = c/n with 



c = rugm/h 2 is the Lieb-Lininger interaction parameter 
and n = N/L is the particle density, to characterize the 
interaction strength between the bosons. For convenience 
we set U = 2m = 1 in the following. 

The Lieb-Lininger model ([TJ is exactly solvable for both 
repulsive (c > 0) and attractive (c < 0) interaction. The 
basic idea is that we first try to find the wave fuction 
in the local region R\ : < x\ < X2 < • • • < xjv < 
L, i.e. assume it takes the form of a superposition of 
plane waves with wave vectors kp j , ip (xi, X2, ■ ■ ■ , xn) = 

Ep, rir . w Apexp V , '•,/••;■ •'•,)• where rj = ±1 repre- 
sent the particle moving to the right or to the left and 
Ap are coefficients to be fixed by the boundary condi- 
tions. Pi, P2, ■ ■ ■ Pn is the permutation of the 1, 2, • • • N, 
is the summation over all the permutations. Un- 
der the open boundary condition, (p (0, x%, ■ ■ ■ , Xn) = 
ip (xi, X2, • • • , L) = 0, the Bethe Ansatz equations are 



exp (i2kjL) = Yl 



i (ki + kj) — c i (ki — kj) + ■ 
i (ki + kj) + c i (ki — kj) — < 



(2) 



with j = 1,2,- • -N. In the region of space under con- 
sideration R : < xi, X2, • • • , Xjv < L, the wave function 
is represented as a summation over all permutations of 
coordinates 

f (xi, • • • ,x N ) = ^2lp(x Pi ,xp 2 , - ■ ■ ,x Pn ) 



x9(x Pl < xp 2 < ■ ■ ■ < x Pn ) . (3) 



with 



6 (xp 1 < ■ ■■ < x Pn ) = 6 (x Pn - x Pn _ 1 ) ■■■9 (xp 2 - x Pt ) 

(4) 

where the step function 8 (x — y) is for x < y, while it is 
1 for x > y. 

The quasimomenta {kj} can be solved from the Bethe 
Ansatz equations @. The energy and the total momen- 
tum of the system are thus E = X^Li k 2 and k — J2jLi kj , 
respectively. Correspondingly the wave function VP is 
known as © , on which all physical quantities can be cal- 
culated. Taking the logarithm of ^ 
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leads us to N — 1 coupled transcendental equations that 
determine both the exact iV-particle ground state and ex- 
cited states of the problem. It is clear that the solution 
of {kj} depends critically on the choices of the quantum 
numbers {Ij}, which is a set of integers, and the initial 
values of {kj}. Here we focus on characteristic properties 
of three typical cases in the current research of ID inter- 
acting cold atoms, all of which are solutions of equation 
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Fig. 2: (Color online) The quasi momentum spectrum of N = 4 
bosons of TG, STG and Bound State branches for interaction 
constant |-y | = 50. 



© with Ij = j. They are: (i) TG Branch: For c > 0, 
the ground state solution of {kj} is a set of real numbers. 
This solution approaches to the TG gas when c — > +oo 
with kj = jn/L. (ii) Bound State Branch: For c < 0, 
the {kj} solution for the ground state is complex which 
is known as the string solution just like the dimer and 
trimer states of bosons with periodic boundary condition 
[26j[27]. For example, the quasimomenta {kj} in the case 
of N = 4 are respectively fcx,2 = ct\ ±*Ai, k^^ = a2±*A2, 
which form pairs with conjugated partners. In the limit of 
strong attractive interaction, the iV-string solution takes 
the form kj = tt/NL + i(N + 1 - 2j)(c/2 + Sj)/L where 
5j is a set of small numbers which exponentially decay 
to zero when c — > — oo and the energy of ground state is 
E/N = k 2 /N 2 L 2 -c 2 (N 2 - 1)/12L 2 . (hi) STG Branch: 
The real solutions {kj} of equation §5§ for c < describe 
one of the highly excited states of the system. In the 
strongly attractive interaction limit c — > — oo, it corre- 
sponds to the STG gas. 

We show the energy spectrum of N = 4 bosons for the 
three branches in the hard wall trap of width L = 1 in 
Fig. 1. As in the case of homogeneous problem of three 
particles with periodic boundary conditions 57J and two 
particles in a harmonic trap |23U28j , the spectrum is quite 
different depending on the sign of the interaction param- 
eter 7 or c (we use — 1/7 in all figures in this paper). 
However, two matching points can be identified for the 
three solutions. In the non-interacting limit the ground 
states for the weakly repulsive (TG Branch) and weakly 
attractive (Bound State Branch) interacting systems join 
smoothly at the left and right end of the axis (7 = 0) 
with E/N — n 2 /L 2 . On the other hand, in the strongly 
interacting limit the TG Branch adiabatically connects to 
the STG Branch at the origin of the axis (— 1/7 = 0) with 
energy E/N = n 2 (N + l)(2N + l)/6L 2 . 

The solutions of quasi momentum {kj} for the three 
cases are presented in the complex plane of k in Figure 
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Fig. 3: (Color online) The transition probability from TG 
branch to STG branch for N = 4 bosons. The system is ini- 
tially prepared in TG gas with 7 = 200 and Lieb-Lininger gas 
with 7 = 1. Inset shows the corresponding transition probabil- 
ity to the BS branch. 



2 for a large interacting parameter I7I = 50. The values 
are real for TG and STG branches, which become closer 
and closer when approaching the strong interacting point 
— 1/7 = 0. Unlike the case of periodic boundary condition 
[20) . they are not any more symmetric about k — 0. By 
assigning the initial value of {kj} a complex number, one 
obtains the string solutions for BS branch with a rather 
small real part and the imaginary parts lie symmetrically 
along the axis Im(k), indicating the conjugacy. 

As shown in Figure 1, the energy changes smoothly from 
TG gas to STG gas when tuning the interaction parame- 
ter — 1/7 across the strongly interacting point from —0 to 
+0. This actually insures the deterministic preparation of 
the STG gas from the TG gas after a sudden change of 
interaction parameter via CIR. However, for the general 
case of an interacting bosonic system, this quench dynam- 
ics is complicated because the resulting state is a linear 
combination of eigenstates of the new Hamiltonian with 
— 1/7 > 0. The transition probability to other states, such 
as the bound states, is non-negligible for weakly interact- 
ing final state. Figure 3 shows the quenching dynamics 
of this interacting boson model, where the initial state is 
prepared in the TG Branch with 7 = 200 or 7 = 1, respec- 
tively. The transition probability to the final state in STG 
branch, characterized as \{^stg I ^tg)| 2 j is plotted for 
varying parameter |7|. The result confirms the metastable 
feature of the STG gas against collapsing into the bound 
state in the strongly attractive interaction regime. On 
the other hand, one has less probability to achieve the 
STG gas if the initial state is weakly repulsive. Inset in 
Figure 3 shows the transition probability to the bound 
state branch {{^bs I ^tg}\ 2 , f° r the same initial states. 
It is more likely that one will arrive at the bound state of 
small I7I if the system is initially prepared in the weakly 
interacting Lieb-Lininger gas. 

Since single-particle density profiles normalized to N 
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Fig. 4: (Color online) The single particle density profiles for the three typical quantum phases of ID interacting cold atoms. Left, 
middle and right panels show the gray scale plots for STG, TG and BS branches respectively, for seven gradually increasing 
interaction parameters and N = 4. Parameters from the top down are 7 = — 0.1, — 1, —2.5, —5, —12.5, —25, —100 for STG, 
7 = 0.1, 1, 2.5, 5, 12.5, 25, 100 for TG, and 7 = -0.1, -0.5, -1, -2, -3, -4, -5 for the BS branches. 



depends crucially on the many-body wave function 'F, one 
might anticipate to more easily recognize the distinction 
between the three quantum phases from their density dis- 
tribution defined as 



p(x) — N I dx-2 ■ ■ ■ I dXN I'F (x, X2 
Jo Jo 



,x N )\ . (6) 



Physically, p (x) represents the probability that the bosons 
occupy the position x, which is accessible experimentally 
by means of the time of flight imaging technique. The 
single-particle density profiles for the three branches so- 
lution are illustrated in gray scale plot in Figure 4 for 
several gradually increasing interaction parameters and 
N = 4. The left panel shows that in the excited state of 
STG branch the density of weakly interacting bosons ex- 
hibits seven peaks due to the slightly attraction between 
atoms, which develops gradually into four in the bottom 
plot when approaching the strong interacting STG regime. 
The density in the ground state of repulsive interacting 
bosons evolves in a quite different way as shown in the 
middle panel, however, merges into an identical four-peak 
structure in the strongly interacting point. This feature 
indicates the fermionization of bosonic atoms for either re- 
pulsive or attractive interaction for large value of |7|. On 
the right panel, atoms in the ground bound state display 
always in a single peak, even for very large interaction 
parameter. We observe again the two matching points at 
7 = and — 1/7 = 0, as can be easily seen from the almost 
identical profiles at the bottom two plots of STG and TG 
branches and those at the top two plots of TG and BS 
branches. 

Another interesting feature of ID exactly solvable model 
is the two-particle correlation function defined as 



g 2 (x) = N(N-l) dx 3 -- 
Jo 

x I* (X,X,X3 • • • , X N )\ . 



dx 



N 



(7) 



which may be used to identify and classify the quantum 
phases of ID Bose gas. It is of particular importance for 
the studies of coherence properties of atom lasers produced 




Fig. 5: (Color online) The two-particle correlation function of 
TG, STG and Bound State branches in unit of g2 (L/2, 7 = 0). 



in ID waveguide, gi (x) denotes the probability that two 
successive measurements will find an atom at the same 
position x. The factor of N(N — 1) in Eq. ([7]) reflects 
there are N(N — 1) ways of choosing the N — 2 integral 
variables out of N, just as we have N ways of choosing 
the N — 1 integral variables out of N in the density dis- 
tribution Eq. ^j. In the Fig. 5, we show the two particle 
correlation function for 4 bosons in the position x = L/2 
by numerically solving the Bethe Ansatz equations and 
doing the multiple integral for wave functions. We em- 
phasize that the method used here is different from the 
approach in [20j . where the correlation function has been 
obtained by numerically solving the integral equation in 
the thermodynamic limit. It can be perceived that gener- 
ally the STG branch is more strongly correlated than the 
TG branch, especially g2(L/2) — 2 for very small value of 
7|. At the strong interacting point, the correlation van- 
ishes for both STG and TG branches, g 2 (L/2) = 0, in 
agreement with the result of periodic boundary condition 
in thermodynamic limit [5D]. On the contrary, the cor- 
relation function for the ground bound state with strong 
attractive interaction is significantly larger than the STG 
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and TG gas, while it decreases rapidly to the result of 
non-interacting system, g 2 (L/2) = 1. 

In summary, we have theoretically studied the ground 
state and excited state of a ID Bose gas with arbitrary 
interaction strength in a hard wall trap. We show in 
the energy phase diagram the three quantum phases TG, 
STG and BS are closely connected to each other in the 
weakly and strongly interacting limits. We also calculate 
the density distribution and correlation function in the 
above cases and conclude that the STG gas is a strongly 
correlated state. 
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